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Abstract 

in 

Dynamics of magnetic bubbles in planar ferromagnets described by the 
Landau-Lifshitz equation with dissipation is analyzed. The pure 0(3) sigma 
C^- | model has static multisoliton solutions, characterized by a number of pa- 

rameters. The parameters describe a finite dimensional manifold. A small 
C$ ' perturbation of energy functional with respect to the sigma model forces soli- 

tons to move. Multisoliton dynamics is effectively reduced to a flow in the 
parameter space. 

a. Bogomol'nyi theories. There exists quite a number of classical field theories with 
the so called Bogomol'nyiB energy bound. The classic example of such a theory, which is 
familiar to both particle physics and condensed matter community, is an 0(3) sigma model 
or a planar ferromagneti described in appropriate units by the energy functional 
rS ■ 

E a [M] = I ( d 2 xd k Md k M , (1) 



2 

where k — 1, 2 runs over planar dimensions and M is a 3-component magnetization vector 
subject to the constraint 

M M = 1 . (2) 

For the energy (0) to be finite the magnetization tends to a constant vector at spatial infinity 
for any time. For definiteness we take (0, 0, 1) as this constant vector. With this boundary 
condition M(t;x) can be viewed as a map from a compactified plane ( equivalent to S 2 ) to 
the S 2 manifold of magnetization defined by the constraint @. The energy functional ([[]) 
is bounded from below 
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E a [M]>4n\N\ . (3) 
N is an integer topological index of the map S 2 — > S 2 , 

N = J d 2 x q(t, x) = — j d 2 x M[diM x d 2 M\ , (4) 

where q(t, x) is a topological charge density. For any winding number N the bound (Q) is sat- 
urated by a static multisoliton configuration characterized by a finite number of parameters 
£. For example a solution with a negative topological index of — n is given by 

w _ .w + jy .W-W 

M B (x;c,ai,bi) = (—— , i- 



■i + iw^i 2 ' i + \w\ 2 ' \ + \w\ 2 ' ' 

- - b n ) 

a's, 6's and c are complex parameters, z = x + iy. Not all parameters are independent, 
some different combinations of parameters have to be identified as they give the same Mb- 
After this identification there are An — 1 real parameters left, they parametrize a (4n — 1)- 
dimensional real manifold M n which will be called a moduli space. It has to be stressed that 
(|) is a time-independent multisoliton solution, its energy does not depend on the choice of 
parameters. 

b. Reiativistic dynamics. The energy functional (0) defines a static version of the sigma 
model. Dynamics can be introduced to the model in a couple of ways. One dynamical 
version is a reiativistic model described by the Lagrangian 

L a [M] = \jd 2 x[ d t Md t M - d k Md k M } . (6) 

A method to study low energy soiiton dynamics in a reiativistic Bogomol'nyi theory has been 
proposed by MantonS and can be briefly summarized as follows. A Bogomol'nyi solution 
like (f|) saturates the lower energy bound @. If solitons are forced by initial conditions 
to move with low velocity, a field configuration at any instant of time remains close to 
the Bogomol'nyi solution and can be approximated by M# [:?;£(£)] with time dependent 
parameters. The approximation is expected to be the better the lower is velocity. Such 
an approximate configuration can be substituted to the Lagrangian (|6]), one obtains after 
integration over plane an effective low energy Lagrangian 

r 2 dM B dM B 
= Jdx ■ (7) 

The field theory ([]) is effectively reduced to a finite dimensional mechanical system. The 
low energy dynamics of solitons is described by a geodesic motion on the moduli space M n 
equipped with the metric tensor g a /3(0- The geodesic approximation has been studied in 
detail for the reiativistic sigma model in §. It has also been explored in other models like 
the original BPS theory of monopolesi or the abelian Higgs model!. 
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c. Landau- Lifshitz equation with dissipation. Another dynamical version of the Bogo- 
mol'nyi theory ([!]), which of main interest for us, is given by 

AT/ 

A M x d t M + T(M) d t M = Pfi[ V 2 M + e^(M) } , (8) 

5M 

where is a projection operator on a subspace orthogonal to M, defined by PA = A — 
M(MA) for any vector A. The equation (§) has to be supplemented by the constraint ([|). 
f (M) is assumed to be a positively definite symmetric matrix for generality eV[M] is a 
small perturbation with respect to the Bogomol'nyi energy (TJ), the total energy functional 
is E = E a + eV. 

Eq. ([D is satisfied for e = by the static multisoliton Bogomol'nyi field (|5]), 

V 2 M B - M B (M B V 2 M B ) = . (9) 

For nonzero £ the solitons are no longer static, they move with velocities proportional to e 

t = 0(e) , (10) 

where the index a numbers the collective coordinates. At the same time magnetization is 
given by the Bogomol'nyi field plus a small deviation 

M(t,x) = M B [x,£{t)} +s m(t,x) + 0(e 2 ) . (11) 

The equation (|8|) does not have Lagrangian formulation in a generic case of T ^ 0, one 
can not proceed along the same lines as in the relativistic case. Instead one has to rely 
on field equations. Eqs. (|IT||TTD define our perturbative expansion in e. It follows from the 
condition (|10D that d t m(t,x) = 0(e). Substitution of Eq.([TT|) to the field equation (D) and 
linearization in e gives 

£ a [ A M B x + r(M fl ) ^ ] + e P, Ib ^(M b ) = 

e [ V 2 rh- M B (M B V 2 m) - (M B V 2 M B )m - M B (mV 2 M B ) ] . (12) 

Similar linearization of the constraint (0) leads to a constraint on m 

M B m = . (13) 

V 2 M B is parallel to M B according to Eq.@. Because of this property and the constraint 
(|I3|) the last term on the RHS of Eq. (|T2|) is zero. 

Eq . fll~2|) is a linear inhomogeneous equation for m. The source term on the LHS of this 
equation depends on the Bogomol'nyi fields only, the RHS can be interpreted as a linear 
operator (dependent on M B [x, £(£)]) acting on the field m(t,x), say, Lm. A projection 

of Eq.(|T2"D on somewhat similar as in i, which is a left zero mode of L, results in a 
solvability condition 

C[\u af i(0 + G Q p(0}-eF(i(0 = 

e J d 2 x — f [ V 2 m - M B (M B V 2 m) - (M B V 2 M B )m } , (14) 
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where the tensors on the LHS are defined by 

,8M B 8M B , 



dM B -» , 8M B 



/ff 



<9£ Q v y 0f0 ' 

MO = -J*x l^Jft$ B ) = ~^J d 2 x V[M B (x,0\ ■ (15) 

If T = 7I with a constant 7, then G a p(£) = jg a p(0 is proportional to the metric tensor (^) 
on the moduli space, as it was discussed qualitatively in 0. Fp(£) can be interpreted as a 
potential force. 

The RHS of Eq.([L4]) is zero. Clearly the second term of the integrand is zero because 
M B ^§r = thanks to the constraint @. After integration by parts the RHS of Eq.flHJ) 
becomes 

e J ** m [ - a -^-(M B V'M B ) ] . (16) 

On the other hand taking a derivative of Eq. (|]) gives 

- ^(Mb^Mb) - Mb^(Mb^Mb) = . (17) 
This equation and the constraint fll3|) imply that the integral (0) is zero. 



To summarize, we have found that a solvability condition for Eqs. ( jl2| , |T3| ) is given by the 
following equation of motion for the collective coordinates 

e[\co a p(0 + G a p(£)]=eFp(Z) . (18) 

Once again the dynamics of a field theory is reduced to a finite dimensional mechanical 
system. 

d. Example. To substantiate the general discussion by a simple example let us consider 
the case of one soliton in an external potential and T = 71 with a constant 7. A general 
form of one soliton solution is given by W = fi/(z — u), where v is a complex position of the 
soliton and the real /1 is soliton's size. Nonvanishing tensor elements are 

uj vV = —ujp u = —2m , 

9vv = 9vv = 2vr . (19) 

g wt is divergent on an infinite plane; it follows from the (5 = \l component of Eq. flT8"D that \x 
is constant provided that is finite. Let the interaction energy with an external potential 
be given by 

V = eJ d 2 x5 i2 \x) q{x,t) , (20) 
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where q is the topological charge density, which is negative in this case. For e > the soliton 
should be attracted by the impurity at the origin. This form of interaction energy appears 
for example in a sigma model for quantum Hall ferromagneti. The potential forces are 

2efi(i/v — fi 2 ) 

- 2eu 2 z/ . , 

F »= F » = , .,2V* • ( 21 ) 



F, 



The equation of motion for v is 



7r(vis + n 2 y 



"-k z [yv + fi 

For 7 7^ the solution of this equation is given by 

R(t) 



Q(t) = 6(0) - (tan B) In 



R(0) 



21n li +3 ^ 2(t) -^ (0)]+ ^^ 2(t) -^ w ' (23) 

where Ae tB = 7r 2 /z 4 (7 — i\)/e, Re l& = vjfi and fi is constant. For A = (purely dissipative 
case) v relaxes to the equilibrium position at v — along a radial line. In general it moves 
towards v = along spiral lines. In the case of Landau-Lifshitz equation or 7 = the soliton 
rotates around the origin along an equipotential circular orbit 

R(t) = R(Q) , 

e(t) = e(0) - wrimp ■ < 24 > 

e. Conclusion. The central result is the equation (|T8D, which gives a prescription how 
to deal with dynamics of solitons in dissipative system close to the Bogomol'nyi limit. The 
equation can be easily generalized to other models because its derivation does not depend 
much on the special properties of the theory (§). 

Generalization to a model with more than one order parameter or field is possible. For 
such a model it would be natural to expect the relaxation times ji for different order pa- 
rameters to be different, 7^ 7^ jj ifi^j. In such a case the tensor G a p is not proportional 
to the metric tensor g a p even if 7's are constants. 

One of applications could be the dynamics of vortices in superconductors at nonzero tem- 
perature. A Bogomol'nyi theory in this case is defined by a Ginzburg-Landau functional for 
a superconductor at a border between type I and type II superconductivity. An appropriate 
small perturbation of the quartic potential!, corresponding to eV, drives the system in the 
direction of weak type II superconductivity. Somewhat similar free energy functional as for 
superconductors describes transition from smectic A to nematic phase of liquid crystals^. 
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